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Abstract. The recent experimental developments in the measurement of hyperfine splittings in the bound states 
of charmonium and bottomonium are presented. Their implications for the hyperfine interactions in the heavy 
quark systems are discussed. 



1 Introduction 



The richness of the spectra of the excited states of atoms 
as well as hadrons lies not only in the principal quantum 
number and angular momentum dependence of the states, 
but in the spin-dependent multiplicities. These arise from 
spin-orbit, tensor, and spin-spin interactions between the 
constituents. Of these three, the most interesting is the hy- 
perfine structure that arises due to the magnetic interac- 
tions between the spins, which causes the splitting between 
spin-singlet (si + S2 - s - 0) and spin-triplet (si + S2 - 
5=1) States. It is the transition between the i and '^o 
states of the hydrogen atom which gives rise to the famous 
21 cm line which is the workhorse of microwave astron- 
omy. 

Hyperfine interaction is equally important in hadron 
spectroscopy. In the quark model, the ground state singlet 
masses of mesons made up of two quarks are simply given 
by 

32;rac , 

M(qxqi) = niiiqi) + mi(q2) + liAWI Si ■ S2 (1) 

ymim2 

where <s,.s.) = -^ for , = 0, =4 for . = 1 
The hyperfine splitting is 

AMhf(nS) = Min^Si) - M(n'So) = ^^l<A«(0)l' (2) 

9m [1712 

It is remarkable how well this textbook prediction works 
with the rather realistic assumption about the strong cou- 
pling constant a'5(M,t/) - 0.6, as{u,d,s) - 0.4, 0-5 (c) = 
0.32, asib) - 0.18, and that |i^(0)P/mim2 is a constant, 
= 33. As shown in Table I, the predicted hyperfine split- 
tings, AM{^S\ 5o) agree very well with their experi- 
mental values. 
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2 Hyperfine Interaction Between Heavy 
Quarks 

The hyperfine splittings of heavy quarkonia require a more 
careful study than their inclusion in Table I might suggest. 

The analogy we have made between atoms and mesons, 
which implies a Coulombic or (oc 1/r) interaction, is not 
quite correct. The qq central potential is known to have 
an additional confinement part which is not well under- 
stood, even though it is generally assumed to be a Lorentz 
scalar, and is parameterized as being proportional to r 

4 as 

Viqq)^-^+Cr (3) 
3 r 

This is the famous Cornell potential, illustrated in Fig. 1. 



2.1 Effect of Quark Confinement 

The Coulombic, 1/r vector part of the Cornell central po- 
tential implies the familiar spin-dependent interactions, the 
spin-orbit, the tensor, and the spin-spin interactions. The 
contribution of the confinement potential to spin-dependence 
does not follow easily. If it is assumed to be Lorentz scalar 
as is normally done, it can not contribute to the hyper- 
fine interaction. However, it may very well have a different 
Lorentz character Only experiments can decide. Further, 
the presence of a confinement interaction poses the addi- 
tional question of how the hyperfine interaction changes 
with the radius of the potential at which the specific states 
reside. This is illustrated in Fig. 1 by vertical lines mark- 
ing the approximate radii corresponding to the |cc) and 
\bb^ bound states. Clearly, the importance of the confin- 
ing potential changes. The interesting questions therefore 
are how the consequent hyperfine splitting, changes with 
the 

1 . principle quantum number, i.e., between 1 S and 2S states 

2. angular momentum, i.e. from L = (S-states) to L 7^ 
(e.g., P-states) 

3. quark masses, e.g., from c-quark to b-qnark states. 
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Table 1. Hyperfine splittings AMi,f = M^S)- M('S) for qq mesons. 



AM,,f{lS) 


(rj, oj) 


(D,D') 




(ric, JW 


('7*,ni5)) 


Eq. 1 (MeV) 


221 


147 


147 


118 


66 


Expt. (MeV) 


234 


142 


144 ± 1 


117 ± 1 


71 ±4 



* With the same parameters the (n,p) and (K,K') splittings are both predicted to be a factor three smaller than observed. This may be 
attributed to the small masses of the Goldstone bosons n and K. 
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Fig. 1. The phenomenological qq Cornell potential. The vertical 
lines show the approximate location of the |cc) charmonium and 
\bb) bottomonium bound states. 

and what the changes imply for the spin-spin hyperfine in- 
teraction between heavy quarks. The best place to address 
these questions is with hidden flavor {q and q of the same 
flavor) mesons, such as the charmonium and bottomo- 
nium mesons. 

2.2 The Experimental Problem 

There is a generic problem in measuring hyperfine split- 
tings, 

AMi,f(nL) = M(n^L) - M(n^L). 

The problem is that while the triplet states are conveniently 
excited in e^e^ annihilation, either directly (e.g., ^Si) or 
via strong El radiative transitions (e.g., ^Pj), the excita- 
tion of singlet states is either forbidden, or possible only 
with weak Ml allowed (n — > «) and forbidden (n — > n') 
transitions. This results in the following situation. 

While the J/iff(l^S i) state of charmonium was discov- 
ered in 1974 |1J and T{l^Si) state of bottomonium was 
discovered in 1977 ||2|. After several false identifications, 
77c(l'5'o) was identified in 1980 |3|, and for more than 
thirty years, the only hyperfine splitting which had been 
measured in a hidden flavor meson was |4| 

AMhf{lS)cc = M(J/i//) - M{i],) ^ 116.6 ±1.0 MeV (4) 



No other singlet states, 

r]'^(2^So)c-c, hc{l^Pi)cc, or J]h{l^So)hh. 

were identified, and none of the important questions about 
the hyperfine interaction which we posed earlier could be 
addressed. 

This has changed in the last few years. 

I want to briefly describe these recent experimental de- 
velopments and their consequence for theory. 

2.3 Hyperfine Splitting in Charmonium Radial 
Excitation: The Search for /7'.(2'5o) 

For charmonium, the hyperfine splitting in this case is 

AMhf(2S),a = Miijj'il^Sx)) - M(i^'^{2'So)) (5) 

The mass of the triplet state radial, if/'{2^S]), is very well 
measured, Miiff'i2^S i)) = 3686.09 ± 0.04 MeV [4J. What 
is required is to identify 77^(2'5o), and to measure its mass 
with precision. The identification of 77', in the radiative de- 
cay 

(A'(2^5i)^r77;(2i5o) (6) 

is diflicult because the transition is a weak Ml, and is pre- 
dicted to have very low energy (Ey ~ 30 - 50 MeV). It 
has never been identified! One has to find other ways of 
populating 77'., for example in photon-photon fusion, or in 
Z?-decays, and to reconstruct it in some of its hadronic de- 
cays. So far only one such decay has been identified. It is 

il',(2S)^ KsKiT. (7) 

The first identification of 77^, came from an unexpected 
source. In 2002, Belle reported 77^ observation in two dif- 
ferent measurements using 42 fb"' of e^e^ annihilation at 
^/s = 10.58 GeV In B decays into K{KsKn) they claimed 
its identification with 56 counts in KsKtt, and M(ri'^) = 
3654 + 10 MeV 111. In double charmonium production, 
Jl^f'^'c they observed a signal with 42^j^ counts, and M(ri'^) = 
3622 + 12 MeV |6|. The need to confirm the discovery of 
77^, and to resolve the discrepancy of 32 MeV between the 
two masses motivated us at CLEO to search for 77^ in two- 
photon fusion with 27 fb"' of data taken by us in the T re- 
gion. We reported observation in KsKn decay with 61;^|5 
events and Mij]'^) - 3642.9 + 3.4 MeV, and consequently 
zfM,,/(25),e = 43.1 + 3.4 f7|. 

Our observation was followed by a similar measure- 
ment by BaBar with 88 fb' of data taken at = 10.58 GeV. 
They observed 112 + 24 events, and reported M{rj'^) - 
3630.8 ± 3.5 MeV fSj. The identification of 77^^ was firmly 
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Fig. 2. Identification of ;7^(2'5o) in KsKn invariant mass, (top) 
from B decays by Belle, (middle) from two photon fusion by 
BaBar, and (bottom) from two photon fusion by CLEO. 



established, albeit with uncomfortably large differences in 
mass. The width of rj'^ remains uncertain within a factor 2 
even to this day. The three measurements are illustrated in 
Fig. 2. Since then both Belle and BaBar have reported more 
mass measurements, and the PDG08 |4| average of all the 
mass measurements is M{ri'^) - 3537 ±4 Me V, which leads 
to the hyperfine splitting 



AMhf{2S)cc = A9 + AM&W 



(8) 



Recall that/JM,,/(15)ce = 116.6 ± 1.0 MeV 0. 

There are numerous pQCD-based predictions for 
AMhji2S)cc, and they range all over the map (and occa- 
sionally even hit 50 MeV). However, it is fair to say that 
nobody expected the 2S hyperfine splitting to be ~ 2.5 
times smaller than the IS hyperfine splitting. A model- 
independent prediction, relating 2S to IS splitting using 
J/t//{lS) and ifr'{2S) masses, and e^e" decay widths, gives 
AMhf{2S)cc = 68 ± 7 MeV, which is also off the mark. 

So far lattice calculations are not of much help. The 
two predictions based on unquenched lattice calculations 
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Fig. 3. Comparing allowed El transitions from tf/'CSi) to 
XcjCPj) states of charmonium with the isospin forbidden tt" tran- 
sition to the singlet P-state hc{'Pi). 



Columbia [9] : zlM,,/(25 )cc = 75 ± 44 MeV (9) 
CP-PACS[10]: zlM,,/(25)ce = 25 =43MeV (10) 

It has been suggested that the smaller than expected 2S hy- 
perfine splitting is a consequence of (/'(25) being very close 
to the \cc} break-up threshold, and continuum mixing low- 
ers its mass, resulting in a reduced difference M{\I/'{2S)) — 
M(ri'^(2S)). However, no definitive numerical predictions 
are available so far. 



2.4 Hyperfine Splitting in Charmonium P-wave: The 
Search for/ie('Pi) 

In this case, we have a very simple, and provocative theo- 
retical expectation, namely 



AMhf(lP) = mCp) - mCp) 







(11) 



This arises from the fact that a non-relativistic reduction of 
the Bethe-Salpeter equation makes the hyperfine interac- 
tion a contact interaction. Since only S-wave states have 
finite wave function at the origin, 



AMhfiL ^ 0) = 0. 



(12) 



We can test this prediction in charmonium by 

- identifying the singlet-P state /ic(l'Pi), and 

- by estimating M(^P), given the masses of the triplet-P 
states 0,1, 2 (^^"0,1,2) 

The experimental identification of /!c( 1 ' P 1 ) is even more 
difficult than that of 77'.. The centroid of the ^Pj states is 
at 3525.30 + 0.04 MeV 14|. If Eq. 11 is true, M(/!,) ^ 
3525 MeV, i.e., ~ 160 MeV below the i/''(25) state from 
which it must be fed. Unfortunately, populating has prob- 
lems. 

- The radiative transition ) — > -yhdl^^) is forbid- 
den by charge conjugation invariance. 

- The only other alternative is to populate in the reac- 
tion ip' — > n'^hc- But that is not easy, because a tt" tran- 
sition (M(7r**) = 139 MeV) has very little phase space. 
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Fig. 4. Illustrating the sequence of decays tf/'(2^Si) 
7T°h,(VPi),h, -» 7?/,(l'5o) used by CLEOto identify A,. 



and further, the reaction is forbidden by strict isospin 
conservation. Nevertheless, this is the only possible 
way of populating he, and we at CLEO had to valiantly 
go for it. 

As with all difficult searches, there is a history of hcCP\) 
searches. In 1982, the Crystal Ball Collaboration searched 
for he in inclusive if/' — > if'hc with 0.9 million if/' . They 
found no evidence of he in the mass region M{hc) = 3440- 
3535 Me V, and established the upper limit n^hc) < 

1.09% (95% CL) 

In 1992, the Fermilab E760 Collaboration, taking ad- 
vantage of the fact that, in contrast to e^e" annihilation, 
he can be directly formed in pp annihilation, searched for 
he m pp ^ he ^ T^J/i//. This time there was plenty of 
phase space for tt", but it was the decay which was isospin 
violating. E760 scanned the region 

^f4pF) = 3522.6 - 
3527.2 MeV with an integrated luminosity - 16 pb"^ 

They reported |12| a statistically significant enhancement 
with ~ 30 counts, which they attributed to he, with 



M(he) = 3526.2 ± 0.15 ± 0.20 MeV. 



(13) 



In 2005, the Fermilab E835 Collaboration, repeated their 
search for he with three times larger luminosity (X. = 48 pb"' 
by combining their data for 1997 and 2000 runs. In the re- 
action pp ^ he ^ TiPJIij/, no evidence for he the region 
in M = 3520 - 3540 MeV was found, in contrast to their 
1992 report. However, they reported lfT3l that in the reac- 
tion, pp ^ he ^ yrje, an enhancement consisting of 13 
counts with significance ~ 3cr, and mass 



M(he) = 3525.8 ± 0.2 ± 0.2 MeV, 



(14) 



In 2005, we at CLEO made the first firm identification 
(significance > 6cr) of he in the reaction 
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Fig. 5. The recoil mass of n° in the decay if/' — » n'^hc. (Top) Full 
and background subtracted spectra for inclusive analysis. (Bot- 
tom) Spectrum of exclusive analysis. 



which is illustrated in Fig. 4. In an analysis of 3.08 mil- 
lion i//' decays he was identified lil4J with N{he) = 178 + 
40 events, 

M(he) = 3524.4 + 0.6 ± 0.4 MeV, (15) 
Slit//' JT^'he) X S2{he yrje) = (4.0 ± 0.8 + 0.7) x 10"^ 

In 2008, we repeated our measurement with 8 times 
larger luminosity, and 24.5 million tfr' (TS). As before, data 
were analyzed in two ways. In the inclusive analysis, the 
photon energy, Ey, was loosely constrained, but the decay 
products of rie were not identified. In the exclusive analy- 
sis, instead of constraining Ey fifteen hadronic decay chan- 
nels of rie were measured. As shown in Fig. 5, a total of 
N{he) = 1282+ 119 events (1146+ 118 from inclusive anal- 
ysis, and 136 + 14 from exclusive analysis) were observed 
with significance > 13cr. Precision results were obtained 

M(he) = 3525.28 ± 0.19 + 0.12 MeV, (16) 
Si X S2 = (4.19 + 0.32 ± 0.45) x 10"^ 



i//' n°he, he 



heCP\) is now firmly estabUshed. 
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If it is assumed that M(^P) is identical to the centroid 
of the triplet-P states, {mCPj)) = [SMixd) + 'iM(xc\) + 
M{Xco)]l9 = 3525.30 ± 0.04 MeV, then the above MQic) 
leads to the hyperfine splitting, 

AMhf{lP)cc = {m(^Pj)) - M('Pi) = 0.08 + 0.22 MeV 

(17) 

But, (M(Vy))^, ^ MCpy. 

The centroid {m{^Pj)^ is a good measure of M(^P) 

only if the spin-orbit splitting between the states 2, ^^"1, 
and ■'f is perturbatively small. It is obviously not so. The 
splitting, MC'Pi) - M(^Pq) = 141.5 MeV, is hardly small. 
Further, the perturbative prediction is that 

mCPi) - mCPo) = ^ [m(V2) - Me Pi)] (18) 
= 113.9 + 0.3 MeV, 

while the experimental value is 

M(Vi)-M(^Po) = 95.9 ±0.4 MeV (19) 

This is a 20% difference! So we are obviously not in the 
perturbative regime. 

This leads to serious questions. 

- What mysterious cancellations are responsible for the 
wrong estimate of Mi^P) giving the expected answer 
that 

AMhf{\P) = 

- Or, is it possible that the expectation is wrong? Is it 
possible that the hyperfine interaction is not entirely a 
contact interaction? 

- Potential model calculations are not of much help be- 
cause they smear the potential at the origin in order to 
be able to do a Schrodinger equation calculation. 

- Can Lattice help? 

2.5 Lattice to the Rescue 

The Coulombic part of the qq interaction is vector, and in a 
non-relativistic reduction of the Bethe-Salpeter equation it 
leads to a contact spin-spin interaction, i.e., it predicts no 
long-range spin-spin interaction. A long-range spin-spin 
interaction can be obtained either by considering a Lorentz 
vector part in the confinement potential, or by considering 
an extension of the one-gluon exchange for the central po- 
tential. 

In a lattice calculation (unfortunately, still quenched), 
Koma and Koma |T6l address the latter possibility. They 
claim clear evidence for deviations from a one-gluon vec- 
tor exchange in the Bethe-Salpeter kernel which leads to a 
5-function spin-spin interaction. They speculate on the in- 
troduction of a pseudo-scalar exchange (a 0"^ glueball?) 
in addition to the one-gluon vector exchange. As uncon- 
ventional as this suggestion is, it is extremely important 
to confirm this. Our measurement of the P-wave hyperfine 
splitting can hopefully shed light on the subject. 



3 Hyperfine Interaction Between Z^-Quarks: 
Thie Searchi for r]hCSo)bb 

The bh bottomonium system is, in principle, the best one to 
study the fundamental aspects of the hyperfine interaction 
between quarks. This is because the ^j-quarks are the heav- 
iest quarks which make hadrons. As a consequence, the 
quarks in bottomonium are far less relativistic ((y/c)^^ ^ 

0.08 compared to {v/c)^g ~ 0.23), and also have smaller 
strong coupling constant (a^^ =0.18 compared to af,~ ^ 

0. 35). Further, as shown in Fig. 1, the bottomonium states 
lie in the qq potential region dominated by the Coulombic 
part, and are therefore least affected by the uncertainties of 
the confinement part of the potential. All these properties 
make perturbative QCD more valid in bottomonium, and 
provide better testing ground for Lattice calculations. 

Unfortunately, until last year we had no knowledge of 
the hyperfine interaction between /7-quarks. The spin-triplet 
T(1^5'i) state of bottomonium was discovered in 1977 [2J, 
but its partner, the spin-singlet 77^(1 '^o) ground state of 
bottomonium, was not identified for thirty years, mainly 
for the same reasons we have mentioned before — the dif- 
ficulty in observing weak Ml radiative transitions. There 
were many pQCD based theoretical predictions which var- 
ied all over the map, v/iihAMhf(lS)b - 35 - 100 MeV, and 
SiY{3S) yrib) = (0.05 - 25) X 10^1 

So, we knew nothing about the hyperfine interaction 
between /7-quarks. 

This has changed now. The ^^(I'S'o) ground state of 
the \bb^ Upsilon family has finally been identified! 

In July 2008, BaBar announced the identification of 
r]b flTl . They analyzed the inclusive photon spectrum of 

r(35) ^ yribilS) (20) 

in their data for 109 million 7(35) (28 fb ' e+e ). BaBar's 
success owed to their very large data set and a clever way 
of reducing the continuum background, a cut on the so- 
called thrust angle, the angle between the signal photon 
and the thrust vector of the rest of the event. BaBar's results 
were; 

Mirib) = 9388.9!^:^ + 2.7 MeV, (21) 
AMhf{lS)b = 71.4!^:^ + 2.7 MeV 
S(r(35) ^ yT]b) = (4.8 ± 0.5 + 1.2) x 10""^ 

Any important discovery requires independent confir- 
mation. At CLEO we had data for only 5.9 million 7(35 ), 

1. e., about 20 times less than BaBar. But we have better 
photon energy resolution, and we have been able to im- 
prove on BaBar's analysis technique. We make three im- 
provements. We make very detailed analysis of the large 
continuum background under the very weak resonance pho- 
ton peaks. We determine photon peak shapes by analyz- 
ing background from peaks in background-free radiative 
Bhabhas and in exclusive decays. And we make a joint 
fit of the full data in three bins of | cos 6t\, covering the full 
range |cos07-| =0-1.0. So, despite our poorer statistics, 
we have succeeded in confirming BaBar's discovery. The 
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results have since then been submitted for pubHcation IfTSl . 
The results agree with those of BaBar. 

Our results are: 



The average of our and BaBar's result for the hyperfine 
splitting is 

{AMhf{lS)b) = M(r{lS)) - Mijih) = 70.6 + 3.5 MeV 



M(/7i) = 9391.8 ± 6.6 ± 2.1 MeV (22) 
AMhfi\S)b = 68.5 + 6.6 + 2.1 MeV 
S(r(35) ^ yr]h) = (7.1 ± 1.8 ± 1.2) x lO^'* 



A recent unquenched lattice calculation predicts (NRQCD 
with u,d,s sea quarks) AMi,f{\S)h - 61 ± 14 MeV. A 
quenched lattice calculation (chiral symmetry and s, c sea 
quarks) predicts AMhf(_lS)b = 70 ± 5 MeV. Thus, as far 



19 International lUPAP Conference on Few-Body Problems in Physics 



Xi,, ISR \ 



3850S09-007a 




0.75 0.85 0.95 

E(Y) (GeV) 

Fig. 7. Illustrating CLEO results for the identification of rjt, in 
a joint fit of data in three bins of the thrust angle, I: | cos drl = 
- 0.3, II: I cos OtI = 0.3 - 0.7, III: | cos Sri = 0.7 - 1.0. 

as the hyperfine splitting for the \bl?j is concerned, lattice 
calculations appear to be on the right track. 

The situation is quite different with the predictions of 
the strength ofthe radiative transition, S( 7(35 ) — > y?7i(15')). 
There are no lattice predictions of transition strengths, so 
far, and pNRQCD predications for forbidden Ml transi- 
tions are orders of magnitude off \2 1 \ . 



4 Summary 

To summarize, we now have well-measured experimen- 
tal results for several hyperfine singlet/triplet splittings in 
heavy quark hadrons: 

|cc) Charmonium: AMi,f(lS) = 116.7 + 1.2 MeV 
AMhf(2S) = 43.2 ± 3.4 MeV 
JMhf(lP) = 0.02 + 0.23 MeV 

\bb) Bottomonium: zfM/,/(15 ) = 70.6 ±3.5 MeV 

In charmonium, we do not have satisfactory understand- 
ing of the variation of hyperfine splitting for the S-wave 
radial states, and for P-wave states 

- For charmonium, we do not have any unquenched lat- 
tice predictions, at present. 

- For bottomonium, lattice predictions are available, and 
they appear to be on the right track. 

- For neither charmonium or bottomonium there are any 
reliable predictions of transitions strength, particular 
for forbidden Ml transitions. 

Much remains to be done. On the experimental front 
it is very important to identify for bottomonium the al- 
lowed Ml transition, T(IS) — > jritilS), and to identify 



the bottomonium singlet P-state, hhCP\). On the theoret- 
ical front one would like to see unquenched lattice calcu- 
lations for charmonium singlets, and, of course, for transi- 
tion strengths. 

I wish to thank the U.S. Department of Energy for sup- 
porting the research reported here. 
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